Charm mixing in the model-independent analysis of correlated D°D° decays 
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I. INTRODUCTION 

Dalitz plot analysis of three-body decays of neutral D 
mesons is a useful tool in various measurements where 
coherent admixtures of D° and D° are observed. This 
technique was initially proposed for the measurement of 
the unitarity triangle angle 7 in B + — > DK + decays 
PJ [2] . Later it was applied to the measurement of charm 
mixing [3J 0] and to the resolution of the quadratic am- 
biguity in the measurement of the angle j3 using a time- 
dependent analysis of the decay B° — > Dir° [5j [6] . Most 
of these measurements are based on the D —> KgTr + ir~ 
decay which offers the best precision among three-body 
D° decays. 

The technique is model-dependent — it depends on 
the complex amplitude of the D° decay which is obtained 
from the D* + — > D°ir + sample using model assumptions. 
The result of the measurement contains therefore model 
uncertainties. In the case of 7 measurement, this uncer- 
tainty (~ 10°) is already comparable to the statistical 
accuracy [3 [8]. 

However, a modification of the analysis is possible that 
allows to perform a completely model-independent mea- 
surement [Tj. It requires the phase space of the three- 
body D decay to be divided into bins. Information about 
the complex phase in each bin can be extracted from the 
quantum-correlated D° decays from ^>(3770) — » DD pro- 
cess. The measurement of the strong phase in bins of the 
D — » Jfg7r + 7r _ phase space was recently performed by 
the CLEO collaboration [9]. This measurement should 
allow to reduce the error of 7 related to the uncertainty 
in the D -> K° s tt + it~ amplitude to 1 - 2°. 

Recently, charm mixing was observed by the Belle and 
BaBar experiments [T01 E] . With degree- level precision, 
the effect of charm mixing can become significant in the 
measurement of 7. It was shown that mixing contributes 
only of second order in the x and y parameters to the 
ADS, GLW and model-dependent Dalitz plot analysis 
methods, and thus can be safely neglected [12]. But the 
effect of mixing on the binned analysis with the phase 
terms extracted from quantum-correlated D° decays is 



of separate interest. 

In this paper, we investigate the impact of the charm 
mixing on the model- independent 7 measurement, and 
show that ignoring the mixing at all stages of the analy- 
sis is safe up to a sub-degree level of precision. We also 
find that in the coherent production of D°D*° system in 
e + e~ collisions, the effect of charm mixing is enhanced 
compared to the case of D°D° production, and propose 
a model-independent method to measure charm mixing 
parameters in time-integrated Dalitz analysis at charm 
factories. The method is sensitive to both mixing pa- 
rameters, x and y, as well as to CP violation parameters 
rev an d otcv- The sensitivity of the proposed method 
can be improved by adding doubly Cabibbo-suppressed 
three-body modes such as D° —> K + ir~ir°. We estimate 
the sensitivity of the proposed method using Monte-Carlo 
(MC) simulation. 



II. MODEL-INDEPENDENT BINNED 
ANALYSIS OF THREE-BODY D° DECAYS 

To introduce the notation we briefly recap the tech- 
nique of model-independent binned Dalitz plot analysis 
of B + — » DK + , D — » Kgir + ir~ decays used to extract 
the angle 7. As usually presented, this does not take 
charm mixing effects into account. 

The amplitude of the B+ -> DK+, D -> K°ir+ir- 
decay can be written as 

A B = A + ree^z+^A, (1) 

where A = A(m 2 Ks7r+ , m% sw - ) = A(m\,m 2 _) is 
the amplitude of the D° — > Kgir + Tr~ decay, A = 
A(mj_, rriL) is the amplitude of the D° —?> Kc;ir + ir~ decay 
(A(m 2 + ,m 2 _) = A(m 2 _,m 2 i _) in the case of CP conserva- 
tion in D decay), rg is the ratio of the absolute values 
of the interfering B+ -> D°K+ and B+ -> D°K+ am- 
plitudes, and 8b is the strong phase difference between 
these amplitudes. The density of the D decay Dalitz plot 
from B + — > DK + decay is given by the absolute value 



2 



squared of the amplitude 

P B = \A B \ 2 =\A + r B e^+^A\ 2 = 

P + r 2 B P + 2\[pl>{x B C + y B S) 



where 



x B = r B cos(5 B + 7) ; y B = r B sin(5 B + 7) . 



(2) 



(3) 



The functions C — C{m\,m 2 _) and S — S(rn 2 + ,m 2 _) 
are the cosine and sine of the strong phase difference 
8d = argA — axgA between the D° —> KgTr + ir~ and 
D" -> K^tt+tt- amplitudes 1 : 

C = cos5 D (m+,m?_) ; S = smS D (m\,m 2 _) . (4) 

The equations for the charge-conjugate mode B~ — > 
DK~ are obtained with the substitution 7 — > —7. Us- 
ing both B charges, one can obtain 7 and 8 B separately. 

In the binned model-independent approach, the Dalitz 
plot is divided into 2J\f bins symmetrically to the ex- 



change m_ -H- mi 



the bin "z" of the Dalitz plot of D from B 



The expected number of events in 

DK+ is 



Ni = hj 



2%/ KiK^i{x B Ci + y B Si) 



(5) 

where Ki is the number of events in the corresponding 
bin of the Dalitz plot of the D meson in a flavor eigen- 
state (obtained using D* ± — > samples) and h B is a 
normalization constant. The bin index "i" ranges from 
—N to N (excluding 0); the exchange m\ O m 2 _ cor- 
responds to the exchange i O — i. The terms and 
Si include information about the cosine and sine of the 
phase difference averaged over the bin region: 



/ \A\\A\ cos S D dV 



Ci = 



"D, 



J \A\UV J \A\HV 



(6) 



Here V represents the Dalitz plot phase space and T>i is 
the bin region over which the integration is performed. 
The terms Si are defined similarly with cosine substituted 
by sine. 

The expected number of events in each Dalitz plot 
bin ( 5 ) is trivially obtained from the probability den- 
sity ( 2 ) by integrating over the bin area, which leads to 
the substitutions P — !• K i: P ->• K-i, C,S -> Cj, 5j. In 
what follows, we only quote the number of events to save 
space. Normalization constants (such as h B in ([5])) are 
also omitted. 

The symmetry under n + f-> 7r~ requires Ci = C_j 
and Si = —S-i. The values of Ci and Si terms can be 
provided by charm-factory experiments operated at the 



1 This paper follows the convention for strong phases in D decay 
amplitudes introduced in Ref. |13| . 



threshold of DD pair production [5]. The wave func- 
tion of the two mesons is antisymmetric, thus the four 
dimensional density of two correlated D 
Dalitz plots is 



K° s n+ii- 



\A C 



= lA^-A^l 2 = 



I 2 2 „./2 

r (m + ,m_,m + , 



P1P2 + P1P2 - 2 V / Pi J P 2 PiP 2 (C 1 C 2 + SiSa) 



(7) 



where the indices "1" and "2" correspond to the two 
decaying D mesons. In the case of a binned analysis, the 
number of events in the region of the (if !j7r + 7r~) 2 phase 
space described by the indices "z" and "j" is 



Mij =KiK-j 



-2sfK~K~K~KZ 3 (C l C ] + SiSj). 



(8) 



Once the values of the terms Ci and Si are known from 
charm-factory data, the system of equations ([5]) contains 
only three free parameters (x B , y B , and hs) for each 
B charge, and can be solved using maximum likelihood 
method to extract the value of 7. 

Note that technically the system ^ can be solved 
without external constraints on Ci and Si for J\f > 2. 
However, due to the small value of r B , there is very little 
sensitivity to the Ci and Si parameters in B + — > DK + 
decays, which results in a reduction in the precision on 7 
that can be obtained [T3] . 



III. CONTRIBUTION OF CHARM MIXING TO 
MODEL-INDEPENDENT 7 MEASUREMENT. 

In the case of CP conservation, the mass eigenstates of 
the neutral D system are given by 



D h2 = -^(D°±D°). 



(9) 



Charm mixing is described by two parameters, xd and 
Ud, which are defined as 



xd 



ui2 — mi 



r 2 -r. 



VD = 



2T 



(10) 



where 7711,2 and Ti^ are the mass and decay widths of the 
mass eigenstates. We use notations x B and yo instead of 
the more common x and y in order not to confuse them 
with the CP-violating parameters x B and y B introduced 
before. The current world average values are: xd = 
(0.59 ± 0.20)%, y D = (0.80 ± 0.13)% M. 
CP violation modifies the expression (19]) to 



£>i, 2 = pD° ± qD C; 



(11) 



where p and q satisfy \p\ 2 + \q\ 2 = 1. CP-violating mixing 
is thus described by two additional parameters rev an d 
o-cv- 



r C ve 



q/p- 



(12) 
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Below we present all the quantities that enter the 
model-independent analysis including the contribution 
of the CP-conserving charm mixing (the corresponding 
quantities are denoted with the prime mark). The full 
formalism including CV violation in mixing is given in 
Appendix [A] 

The number of events in the each bin of the flavor- 
tagged D° — > KgTr + Tr~ Dalitz plot after integration over 
decay time is 

K[ =Ki+ y/KiK-iivuCi + x D Si) + O (x 2 D ,y 2 D ) . (13) 

Similarly, one can obtain the number of events for the 
D — s- Kgir + TT~ decay from B + — s- DK + : 



N[ =N { + y/KiK-i {y D Ci + x D S z ) 
+r|V 'KiK-i {y D Ci - x D S t ) 
+Ki (x B y D - VbXd) 
+K-i {x b vd + Vbx d ) + O (ijj, y 2 D ) 



(14) 



The amplitude for the correlated D°D° decay, assum- 
ing that the particle denoted with the index "1" decayed 
first, equals 

A corr (t, m\ , m 2 _ , m' 2 , m' 2 ) = A x {0)A 2 (t) -A 1 (0)A 2 (t). 

(15) 

Before the particle "1" decays, the amplitude stays an- 
tisymmetric and the mixing does not affect it. We can 
therefore assume that the particle "1" decays at the time 
t = and count time from the moment of its decay. 

After the integration over time and taking the abso- 
lute value squared of the amplitude (fTHI), we obtain the 



number of events My in the bin u ij" of the phase space 
for the pair of D mesons (still assuming that the particle 
"1" denoted here with the index "z" decayed first): 



(16) 



-2^/KiK-iKjK-jiCiCj + SiSj) 
-Kj^KiK-i(y D Ci - x D Si) 
-K-j \J KiK_i(y D Ci + x D Si) 
+Ki^f Kj KZj (y D C 3 - x D Sj) 
+K- i ^K j K- j (y D C j + x D Sj) 
+0{x 2 D ,y 2 D ). 



However, measurement of the decay time in an experi- 
ment with symmetric beams is difficult. Therefore one 
has to average over the decay order, which leads to the 
cancellation of all terms linear in mixing parameters: 



= Mij 



0(x 2 D ,y 2 D ). 



(17) 



The real analysis performed at CLEO uses the values 
M!-j (which, as we have seen, are unaffected by mixing 
at first order), and values of K{ obtained from correlated 
DD decays where one of the D mesons serves as a flavor 
tag. The Ki values extracted this way are also unaffected 
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by mixing. Therefore, the values of Cj and Si extracted 
in this analysis contain no linear mixing contribution. 

As far as processes observed at B factories are con- 
cerned, both the D — > if <j7t + 7t - decay from B + — > DK + 
and the flavor-tagged D — > A^7r + 7r _ decay contain mix- 
ing contributions, and therefore the observable numbers 
of events in the Dalitz plot bins are N[ and K' i: respec- 
tively. Clearly, if one uses the values Ki, obtained from 
a charm factory, in the fit to obtain 7 from described 
by Eq. |14| the resulting value contains contribution in 
first order in xu, yo- If the values K[ are used, Eq. 
can be rewritten as 

N[ = K[W B Ki l+ 2^K[K'^BC[+y B S^+0{x 2 D ,y 2 D ) , 

(18) 

i. e. it has the same form as Eq. [5j without mixing 
contributions up to second-order in xd, Ud, but the phase 
terms S| contain first-order mixing corrections: 

K '■ + K' K' — K' 

CI = a + ^=(1 - C 2 )y D + -^=^CAx D , 

K' - K'_ K' + K'_ 

S ' 1 = Sl ~ IjEJK^} 1 ~ S * )XD ~ JKTK^ CiSiVD - 

(19) 

Thus, if the terms Ci, Si are left as free parameters in 
the fit to B decay data, the mixing correction is only of 
second order (the effective parameters C-, S"- are mea- 
sured in this case [H]), but if these terms are obtained 
from correlated D decays, first order mixing corrections 
to 7 appear. However, these corrections are additionally 
suppressed by factor r B ~ 0.1, and the residual contri- 
bution of charm mixing to x B and y B is at a percent 
level. 

A quantitative estimate of the effect was performed us- 
ing the procedure described in Appendix [B] Three anal- 
ysis strategies have been considered: 

1. Using Ki (unaffected by mixing) from the coherent 
DD production, 



2. Using K[ (Eq. 13) measured in D*+ -> D n + , D 



Kgir + n decays, 

3. Using Ki and applying a linear correction for the 
mixing contribution according to Eq. 
that xd, yo are known). 
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(assuming 



The effect of mixing on the fitted value of 7 depends on 
the value = arctan(y £>/;££>), the ratio \Jx 2 D + y D /r B , 
and the 8 B and 7 values. In our study, we use 
\J x 2 D + y D l t b =0.1 (all biases are proportional to this 
quantity) and scan over the other parameters. The re- 
sults are shown in Table [i] Clearly, if the K[ values from 
D* + — > D°ir + , D — > Kgir + Tr~ decays are used in the 7 
fit, the mixing contribution can be neglected. 

For clarity, we would like to compare our result with 
the conclusions of previous papers that have addressed 
the impact of charm mixing on 7 measurements [121 116) . 



4 



TABLE I: Estimates of the charm mixing effect on the 7 value 
measured using the model-independent time-integrated Dalitz 
plot analysis method, for three different analysis strategies. 
The maximum 7 biases for varying old = &rctan(yD /%d), Sb 
and 7, as well as the values of these parameters at which 
the maximum is reached, are shown. The estimation assumes 
y/x 2 D + y 2 D /r B =0.1. 

Strategy A7 max Qmax OB, max 'YSjinax 

1. Using Ki 2.9° 184° 85° 87° 

2. Using K'i -0.2° 97° 2° 90° 

3. Linear correction 0.07° 324° 72° 73° 



TABLE II: Statistical sensitivity to the mixing and CP viola- 
tion parameters for the time-dependent Dalitz plot analysis. 
Two strategies are considered: (i) Ki fixed from charm factory 
data and (ii) Ki taken as free parameters. 



Parameter 


Precision 


Ki fixed 


Ki floated 


xo (10- 4 ) 


17 


22 


VD (10- 4 ) 


13 


16 


rev (10~ 2 ) 


9 


9 


oicv (°) 


5 


5 



Ref. [16] considers the case when the D amplitude does 
not contain the mixing contribution (or is corrected for 
it), but the B decay data are uncorrected for mixing. 
This corresponds to the analysis strategy 1 in our MC 
study. The systematic bias to 7 in that case is linear 
in xd and yu and can be numerically large. The treat- 
ment in [12] corresponds to the case when the mixing is 
neglected in both the flavor-tagged D and B data; the 
systematic bias in 7 is second order in X£> and yo in that 
case. In the context of the model-independent binned 
Dalitz plot analysis, the conclusions of [T2] can only be 
applied if the phase terms Cj and S% are left as free pa- 
rameters in the fit to B data. The analysis procedure 
considered here is an intermediate case: part of infor- 
mation about the D — > K i §it + 'k~ amplitude (namely, Ki 
terms) is extracted from the flavor-tagged data uncor- 
rected for mixing, while the Cj and Si terms from the 
quantum-correlated DD decays contain no mixing con- 
tribution. This results in a bias linear in xo,yD (but 
numerically small due to additional rs suppression). 

IV. TIME-DEPENDENT CHARM MIXING 
MEASUREMENT 

The charm sector is the only place where contributions 
to CP violation from down-type quarks in the mixing di- 
agram can be explored. While values of charm mixing 
parameters are not easy to predict in the SM, CP vio- 
lation in mixing is expected to be very small. However, 
there is a range of SM extensions which allow sizable CP 
violation effects [TTJ [TB]. Because of that, precise mea- 
surements of mixing as well as CP violation parameters 
are essential. 

The most accurate measurements of the charm mixing 
parameters to date have been performed by £?-factories 
using time-dependent methods |101 lllj . For example, 
by observing the wrong-sign decay D° K + tt~ [TTj . 
BaBar determines Rjj, the ratio of doubly Cabibbo- 
suppressed (DCS) to Cabibbo-favored (CF) decay rates, 
and the mixing parameters x'p and y' D , where x' D = 
x D cos 5 Krr + y D sin 5 K „ , y' D = ~ X D sin 5 K tt + y D cos 6k-k 
and Sktt is the strong phase between the DCS and CF 
amplitudes. Since only y' D contributes linearly to the 



time dependence of the decay rate, and the strong phase 
6 Kir is close to zero [19] . this kind of measurement is 
practically insensitive to the parameter xp. 

The Belle collaboration has performed a measurement 
of mixing parameters using time-dependent Dalitz plot 
analysis of the D° — > Kgir + ir~ decay [3]. A similar anal- 
ysis has been performed by the BaBar collaboration us- 
ing in addition D° -> K%K+ K~ [4]. Here the Dalitz plot 
distribution depends linearly on both mixing parameters. 
However, model assumptions in the decay amplitude de- 
scription are necessary in this analysis, which unavoid- 
ably results in significant model uncertainties. 

The model-independent binned Dalitz plot analysis 
may be extended to time-dependent measurements of 
charm mixing parameters and CP violation in charm mix- 
ing. The time-dependent number of events in the bin "£" 
of D° —> K s tt + tt~ decay Dalitz plot is 

^ (t) = h D e- rt [K t + \J K i K_ l {y D C % + x D Si) Tt+ 
0({x D +y D f(Ttf) . 

(20) 

Using parameters Ci and Si determined from indepen- 
dent measurements at a charm factory, we can eliminate 
completely the model uncertainty in xd and yo extrac- 
tion. 

We estimate the statistical sensitivity of the time- 
dependent fit to the mixing and the CP violation pa- 
rameters. Note that the values of Ki can be measured 
independently with high precision at charm factory ex- 
periments. However, in that case one has to deal with 
systematic uncertainties due to detector efficiency and 
other effects. Depending on the magnitude of these un- 
certainties, one can consider a fit with Ki as free param- 
eters, where the statistical uncertainty increases, but the 
systematic effects are minimized. 

Table |H] shows the results for both strategies. The toy 
MC simulation uses 10 6 flavor-tagged D mesons from D* 
decays corresponding to the samples available at the B- 
factories; the fit with fixed Ki uses in addition 2 • 10 6 
flavor-tagged D decays in C = — 1 state that can be ob- 
tained at a future charm factory experiment. The sim- 
ulation uses the amplitude of D —> KgTr + ir~ decay de- 
termined by Belle [7], and the binning of the Dalitz plot 
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with 8 bins defined by the uniform division of the strong 
phase difference So [13]. The values of the phase terms 
Ci, Si are calculated from the D — > Kgir + n~ amplitude; 
the contribution of their statistical precision in the final 
result is negligible with the current experimental data. 
Our estimates do not account for uncertainties in the 
time measurement and background effects. 



V. TIME-INTEGRATED CHARM MIXING 
MEASUREMENT 



TABLE III: Statistical sensitivity to the mixing and CV vio- 
lation parameters for the time-integrated Dalitz plot analysis 
using a 10 3 fb _1 e + e~ — > V>(4040) data sample. 



Parameter 


Precision 


Coherent only 


Incoherent only 


Both 


X D (10- 4 ) 


12.5 


18.4 


11.8 


yo (i(r 4 ) 


8.7 


12.9 


8.5 


rev (KT 2 ) 


5.4 


5.2 


3.8 


acv (°) 


3.5 


3.5 


2.5 



The cancellation of terms linear in ifl and yu for 
the correlated decays occurs only in the case of the an- 
tisymmetric wave function of two D mesons produced 
with the charge conjugation quantum number C = — 1. 
It is possible, however, to produce a pair of D mesons 
with both quantum numbers C = ±1 in the process 
e+e- 0(4040) -> D°D*° [20]. 2 The e+e" ->■ cc 
cross-section at the -0(4040) mass is dominated by this 
process and is comparable to that at ■0(3770) [21]. De- 
pending on the D*° final state (D°ir° or D ^), the DD 
pair is either C = — 1 or C = +1, respectively. In the case 
of C = +1, the amplitude of the decay is 



for C = +1, and 



A X A 2 + A X A 2 



(21) 



Unlike Eq. [17] now we allow two D mesons to de- 
cay into different final states, therefore we denote the 
numbers of events in bins of flavor-specific decays as 
ki and Ki, and the corresponding phase terms as Ci,Si 
and C i: Si. These can be either two-body (in that case 
only one bin makes sense, and the index "i" only takes 
the value ±1) or multibody final states. For non-self- 
conjugate final states, the index takes positive values for 
D° decay and negative values for D° decay. The terms 
related to the strong phase difference are defined simi- 
larly to Ci and Si in Eq. [6] The number of events in 
phase space bins with the contribution of mixing taken 
into account is 



— k^ j ~\~ k ^ "^^j7 



+2 i/fcj k-i Kj K^j (aCj + s,Sj) 
+2K j ^k i k^ l (y D c i - x D Si) 
+2K- j ^kikZ l {y D c i + x D Si) 
+2ki y/ KjK_j (y D C 3 - x D Sj ) 
+2k-i^/K j K- j (y D C j + x D Sj) 
+0(x 2 D ,y 2 D ) 



(22) 



2 The optimal energy for this kind of experiment is 4.01 GeV, for 
which the DD* cross section is maximized, and D*D* produc- 
tion is below threshold. 



M~/ =hK-j + k. 



~2\/kik—iKjK—j (ciCj + SiSj) 
-0(x 2 D ,y 2 D ) 



(23) 



for C = -1. 

Important special cases for the two-body decays are: 



1. Flavor-specific decay (e. c, 
1, K-i = 0, 

2. D a -> K-7T+ decay: K x 
cos <5d , S\ = sindjj. 



D { > 



K-e+v e ): K x 



1, K_t 



Ktt> 



The effect of mixing is linear in the case of C = +1, 
which allows to measure the parameters xd and yu in the 
combined fit of the decays of both charge parities. The 
simplest strategy involves flavor-tagged D — > Kgir + ir~ 
decays produced in both C = +1 and C = —1 states. 
Decay of the C = — 1 state yields the numbers of events 
in D — > Kgir + ir~ phase space bins unaffected by mixing, 
while the C = +1 state contains linear mixing contribu- 
tion: 



Mf = + 2^/K l K^(y D C l - x D S % ) + 0{x 2 D ,y 2 D ). 

(24) 

Since there are bins with \Ci \ ~ 1 and bins with \Si\ <~ 1, 
the sensitivity to both xjj and yjy should be of the same 
order. Note that the mixing term is twice larger with 
respect to that in D* + — s- D°ir + , D —> Kg7r + ir~ decays 
(Eq. _ 



13) 



We have made a quantitative estimate of the sensi- 
tivity to mixing parameters using MC simulation. The 
results of this study for a data sample equivalent to the 
integrated luminosity of 10 3 fb _1 are shown in Table III 



This corresponds to one standard year of data taking 
with a peak luminosity of 10 35 cm _2 s _1 . Projects of 
cr-factories with such luminosity are being considered 
|22) . The numbers of events are estimated using the 
double-tagged event efficiency obtained by CLEO in the 
reconstruction of ^(3770) decays [H] (see Table [V] in Ap- 
pendix |b| . 

Our study involves only the flavor-tagged and double- 
tagged D —> KgTr + ir~ decays in the coherent C = — 1 
state (to extract Ki and Ci,Si, respectively), as well as 
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the flavor-tagged D 



decays in the C 



state and the incoherent flavor-tagged D — > Kgir + ir~ de- 
cays for the final extraction of mixing parameters. The 
incoherent events can be obtained at a cha rm factory 
from ^(4040) -> D ± D*T decays. Table [Till shows that 
adding the incoherent decays to the coherent sample 
does not improve the mixing parameters precision sig- 
nificantly. It is explained by the fact that the precision 
is determined by the flavor-tagged D — > Kgir + ir~ decays 
in C = —1 state. However, the inclusion of incoherent de- 
cays improves the precision of CP violation parameters. 

A more complicated analysis can also include the 
double-tagged D — > K s ir + Tt~ decays in C = +1 state 
which also include the mixing parameters at the first or- 
der. Our estimation of the mixing parameters precision 
is based on the assumption that systematic errors can be 
essentially suppressed since both C = +1 and C = — 1 
decays have similar kinematics and may be detected si- 
multaneously in the experiment. 



Equations ( 22 ) and ( 23 ) imply the absence of direct CV 



violation in D decays. It is, however, possible to take its 
effect into account. This requires doubling the number of 
parameters Ki, Cj, Si (by treating decays of D° and D° 
separately for self-conjugate final states like Kgn + TT~, 
or by decoupling the numbers of events in CP-conjugate 
modes for final states like K~ir + ir°). The number of 
equations used to constrain the phase terms Ci,Si from 
the double-tagged D — > KgW + 7T~ decays will stay the 
same, but for a reasonably large number of bins the sys- 
tem of equations will remain solvable. Consequently, the 
method can be used to distinguish between direct CP 
violation and CP violation in mixing. 

CV violation in mixing of neutral kaons can mimic CV 
violation in charm if the final state Kgir + ir~ or other 
states including Kg mesons are used in the analysis. CP- 
violating terms in the amplitudes Ai% are of the order 
eA 2 , where e ~ 2.2 x 1(T 3 23J is the magnitude of CV 
violation in kaon mixing, and A 2 = sin 2 8q — 0.23 2 gives 



the relative difference between the D° 



and 



D — > K^tt tt~ amplitudes. If not accounted for, CV vi- 
olation in the kaon sector will introduce fake CV violation 
in charm of the order r C p ~ e\ 2 /(x 2 D +y 2 D ) ~ 1%. How- 
ever, once the amplitude of the D° — > K\tx + -k~ decay is 
measured (this can be done at a charm factory exploit- 
ing kinematic reconstruction of the missing K\ at DD 
threshold), this effect can be corrected for in the charm 
mixing measurement. 



VI. EXTENSION TO OTHER D° DECAY 
MODES 

The precision of the charm mixing parameters mea- 
surements can be improved by adding other three and 
four-body hadronic final states such as K~tt + it and 
K~TT + ir~Tr + in addition to D — > K s tt + tt~ . In these cases 
the phase terms Ci, Si can be defined in the same way 




m 2 Kji+ [GeV 2 /c 4 ] 



FIG. 1: Uniform phase binning obtained from the CF and 
DCS D° -» K^tv+tt decays amplitude models. 



The advantage of using flavor-tagged Cabibbo-favored 
and doubly Cabibbo-suppressed D° — > K^^ir three- 
body decays for the mixing measurement is that the rel- 
ative size of the interference term can be made much 
larger than in the case of D — > Kgir + ir~ decays. For ex- 
ample, if we observe decays where one of the D mesons 
is tagged by its semileptonic decay to be D°, and the 
other one is reconstructed in the K~tt + tt° state ("wrong- 
sign" tag), the first term in Eq. 24 is of the order 
Rktth ~ 0.06 x 0.06, while the interference term is of 
the order \ZRk-k-k{xdtVd) ~ 0.06 x 0.01. While the sta- 
tistical sensitivity to the mixing parameters should be 
similar compared to the previous analysis, the system- 
atic error can be reduced since the relative magnitude of 
the interference term that contains the mixing parame- 
ters is larger, and thus is less sensitive to background and 
efficiency uncertainties. 

Since the statistical sensitivity depends on the struc- 
ture of the decay amplitude, we perform a MC study 
for the D° — > K^tt 7r° final state. The description of 
CF and DCS D° — > K^^tt decay amplitudes is based 
on quasi two-body models obtained by CLEO [24] and 
BaBar gS]. 

As in the case of the K'Itt + 'k~ decay, we use a uniform 
phase binning |13) : Fig. nl shows the binning obtained 
from the CF and DCS models used in our study. We use 
a sample corresponding to an integrated luminosity of 10 3 



flr 



as for the D — > KgTt + ir decay). The results of our 



as for D — s- Kzir 



calculation of the statistical uncertainties of the mixing 
parameters as functions of the average strong phase dif- 
ference between CF and DCS decay amplitudes 5^ n7T 
are shown in Fig. [2j Apparently, there is a symmetry: 

<t*d ( S %™°) = oyn (&%™° ±k/2). The strong phase 



difference obtained by CLEO equals 227^ \ 7 ° . Our esti- 
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A model-independent method to measure charm mixing 
parameters from time-integrated analysis at a charm fac- 
tory is proposed. Using a MC simulation we find that 
the sensitivity to both mixing parameters xn. and yo is 
of order 10" 3 for one year of data taking with luminos- 
ity 10 35 cm^s- 1 at the peak of ^(4040) resonance. The 
method does not rely on absolute branching fraction mea- 
surements, and therefore does not contain uncertainties 
related to measurements of absolute efficiency or values 
of the branching fractions of the decays involved. Due 
to strong phase variations over the phase space of the 
D — » Kgir + ir~ decay, nearly equal sensitivity to both 
xd and yn is obtained. The proposed method can also 
be used to measure the CV violation parameters of charm 
mixing. 



FIG. 2: Statistical precision for id (triangles) and yD 
(squares) as functions of the average strong phase difference 
between CF and DCS D° -> K^-k+tt decay amplitudes. 



TABLE IV: Statistical sensitivity to the mixing and CV vi- 
olation parameters using D — > Ktytt decays. Results of toy 
MC simulation for the integrated luminosity of 10 3 fb _1 and 
& kW = 227 o^ 



Parameter 


Precision 


Coherent only 


Incoherent only 


Both 


x D (lir 4 ) 


6.2 


8.7 


6.0 


yn (icr 4 ) 


6.1 


8.9 


6.1 


rev (10~ 2 ) 


2.6 


2.2 


1.7 


acv (°) 


2.4 


2.4 


1.7 



mates of the mixing parameters' precision for this value 
are shown in Table IIV1 
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Appendix A: Formalism for coherent D° decays with 
charm mixing 

We denote the amplitude of D° (D°) decay to the final 
state / as A (A). In the case of decay to three spin- 
zero particles, an amplitude A depends on two kinematic 
parameters (here we use two squared invariant masses 
m\ 2 and m\ 3 of the decay products) . The density of the 
D° decay Dalitz plot in the absence of charm mixing is 



VII. CONCLUSION 

We have investigated the impact of charm mixing on 
all stages of the binned model-independent measurement 
of the CKM phase 7 from B+ -> DK+ , D -> K s n + ir- 
decays. We show that ignoring the mixing at all stages of 
the analysis is safe at the current statistical accuracy: the 
bias on 7 with the currently measured mixing parameters 
is of order 0.2°. 

A model-independent approach to perform a time- 
dependent Dalitz plot analysis of the D — > K^K+ir^ de- 
cay to extract charm mixing parameters was considered. 
The approach is promising for high-statistics analyses at 
the LHCb experiment and at SuperB factories. 

We find that in the coherent production of the D°D*° 
system in e + e~ collisions, the effect of charm mixing can 
be enhanced if the D*° is reconstructed in the D°7 decay, 
which requires the D°D° system to be in a C — +1 state. 



P (mi 2 , m 2 13 ) = \ A (m? 2) mf 3 ) | 
Charm mixing modifies this expression as 



P' (m\ 2 ,m\ 3 ,t 



x(t)A + i-a{t)A 
P 



(Al) 



(A2) 



where t is the time between the production and the decay. 
Here we omit the dependence of A and A on the Dalitz 
plot variables. The time-dependent terms are given by 

Tt 

x(t)+io-(t)=exp—(-l + x D -iy D ). (A3) 



After integration over time we get 
P' (m 12 ,m 13 ) =a P + a 1 r^. v P+ 



(A4) 
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where 



r C ve 



J 

P 



S ± = sm(A5 D ±a C v), 



(A5) 



cos (AS D ± a C v), 



ASjj = (5d (Tni2: m i3} — $d ( TO i2) m i3)) i s a difference 
of strong phases for the decays D° —> f and D° —> /, 
and 



a 
ai 

a 3 



1 



1 

1 + 



-x 2 D +yl) +0({x D +y D ) 3 ) 



l-yl l + x 2 D 



(A6) 



= o (*d + Vd) +0((x D + y D ) 3 ) 



Vd 



VD + 0({x D + y D f) , 



T^X = ^ + 0((.T D + y D ) 3 ) 



substitutions p o q and P O P: 



P (ml 2 ,ml 3 ) =a P + a 1 r c £P- 



r^VPP (C-a 2 + S-a 3 ) 



(A7) 



Now we consider the decay of the coherent D°D° pair 
in a C = — 1 or C = +1 state. Assuming the particle 
denoted with the index "1" decayed first, the Dalitz plot 
density accounting charm mixing effects is given by the 
expression 



The expressions for D° decays can be obtained after the 



where 



P corr (( TO 12)i , ("^2)2 > ( TO 13) 2 ) = 6 



P X P 2 + PiP 2 + 2CdP 1 P 1 P 2 P 2 {C X C 2 + 5i5 2 ) 



CV J 



P1P2 + r 2 cv P 1 P 2 + 2C\jP x P x P 2 P 2 (CtC^ - S+Sl 



-b C 3 



P 2 P 2 C+ (r cv Pi + r^Pi) +C^/P 1 P 1 C+ (r cv P 2 + r^P, 



P 2 P 2 S+ (r CP P 1 - r^Pi) +Cy/P 1 P l S+ (r CP P 2 - r^P 



(A8) 



6f 



l+Cyl 



l + Cy 2 D 



(1 + C) yp r 

(l + C)x D 
(1+^) 2 ' 



Cx D 



(1 + ^) 2 



1 



Cx D 



(1 + ^) 2 
(1+C)a 2 , 

(1 + C)a 3 , 



a o + — ^— (-x D + y D ) 



{C + 2) ai , 



(A9) 



and C = ±1 for the symmetric and antisymmetric states. Note that in the case C = —1 the interference term 
vanishes [57] , while for C = +1 it is doubled compared to the incoherent case. 

The expression for the D —> KgTr + ir~ decay density in the B — > DK ± process accounting for the charm mixing 
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contribution is 



,± 



P' B ± {m ± ,m T ) = a 
+ a\ 
+ a 2 
+ a 3 

where x B = r B cos (8 B ± 7 ± «cp), J/s = r B sin (£ B ± 7 ± «cp) 



P + r 2 B P + 2VPP {x B C + y B S) 



,%P + r^r B P + 2VPP(x R : ± C ± -y± ± S ± 



X B ± [ r CvP 



r 



'PP r, 



±i 
CP 



y B i 

r CV r B, 



(A10) 



( r cv ~ r c-p r s) <5 



TABLE V: Numbers of events used for mixing parameters 
measurement. The numbers for K~tt + tt° include both CF 
and DCS decays. 



Type 


Number of events (10 6 ) 




K° S -K + -R- 




Incoherent flavor tagged 


6 


30 


Coherent C = — 1 flavor tagged 


2.1 


10.5 


Coherent C — +1 flavor tagged 


1.4 


7 


Coherent C = — 1 


0.6 
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ment is studied, the values of Xi are fixed (which cor- 
responds to infinite statistics), while in the study of the 
statistical sensitivity of the mixing measurement in Sec- 
tions [V] and VI we sample the Xi values according to the 
Poisson distribution. 

The expected numbers of events (Xi) include the free 
parameters of interest (x B and y B in Section III or mixing 
parameters in Sections [v] and VI ) and are expressed via 



the numbers of events in flavor-specific D° state iQ and 
the phase terms Ci, Si. The values of Ki, d and Si are 
calculated from the D — > K^-k+ii^ and D° 
decay amplitudes [71 [Ml US] as 



Appendix B: Numerical calculation procedure 



The numerical results presented in Sections III V]and 
VI are obtained using a procedure described below. 

In the model-independent approach, the phase space of 
the decay is divided into bins and we deal with numbers 
of events in each bin. The maximum likelihood method 
is used to obtain the parameters of interest from these 
numbers. The likelihood function is defined as 



-21og£=-2^1ogP(X i ,<X i }), 



(Bl) 



where P(X, (X)) is the Poisson probability function to 
observe X events with an expected number of events (X) . 

The observed numbers of events Xi are obtained using 
the exact formulas with mixing (see Appendix |Aj. When 
the systematic effect of charm mixing on the 7 measure- 



a = 



s, 



Ki = pDdm 2 + drn 2 _, 



n 



J v% y/PDPD C dm\dm 2 _ 
J v p B dm 2 ^dm 2 _ J v p^dm^dm 2 . 



J Vz VpdPd S dm 2 + dm 2 _ 
J v ppdm+dm?, podrn+dm 2 . 



(B2) 
(B3) 

(B4) 



The free parameters are extracted using MINUIT to min- 
imize the likelihood function (Bl I. 



Table [V] shows the numbers of events that we use to 
obtain the results in Sections IVl and IVT1 To estimate the 
numbers we use the detection efficiency determined by 
CLEO [9] and the statistics corresponding to an inte- 
grated luminosity of 10 3 fb _1 at the -0(4040) resonance. 
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